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1. INTRODUCTION 

The issue of color confinement in quantum chromodynamics [QCD], widely accepted as the theory of the strong 
interaction, is a longstanding problem both conceptually and quantitatively. One necessary prerequisite for color 
confinement is the conservation of color charge: a 'leaky' system clearly cannot be confining. Within the framework 
of functional field-theoretic methods in nonabelian theories (i.e., QCD), expressed in terms of local Green's functions, 
the Slavnov-Taylor identities [1, 2] are the consequence of this charge conservation. 

The Slavnov-Taylor identities come in various guises and have been applied to many problems, in particular to the 
nonperturbative study of the Dyson-Schwinger equations. For example, in quantum electrodynamics [QED] where the 
Slavnov-Taylor identities reduce to the Ward identity, their study led to the Ball-Chiu vertex [3] with its subsequent 
improvement, the Curtis-Pennington vertex [4], which have been used to study dynamical mass generation (see for 
example Refs. [5, 6]). In axial gauge Yang-Mills theory, the Slavnov-Taylor identities [7] have been used to study the 
gluon propagator Dyson-Schwinger equation, with the inference that the propagator diverges as l/q^ in the infrared 
(one potential signal for confinement) [8]. This result is however, not without ambiguity [9]. In Landau gauge Yang- 
Mills theory, the Mandelstam approximation [10] (see also Ref. [11]), where ghost contributions to the three-gluon 
vertex Slavnov-Taylor identity are neglected, one can also obtain an infrared enhanced gluon propagator. The above 
examples can be colloquially referred to as Schwinger-Dyson studies (these are early studies and a quirk of history 
has subsequently resulted in the reordering of the names) and they share two common features: the utilization of the 
Slavnov-Taylor identities and either require no, or assume no ghost contributions. 

In attempting to further the Mandelstam approximation to Landau gauge Yang-Mills theory by including the ghost 
contributions, a startling new phenomenon emerged — infrared ghost dominance [12] (see also [6, 13] for contemporary 
reviews). The starting point was an approximation to the Slavnov-Taylor identity for the ghost-gluon vertex. This 
approximation was shown to be inconsistent with perturbation theory [14], but this turned out to be unimportant: 
using general arguments one can show that the tree-level part of the ghost-gluon vertex is the important ingredient 
[15]. Indeed, there are indications that the dressing of the ghost-gluon vertex function is only slight [16]. The results 
of this infrared ghost dominance can be summarized as follows: the gluon propagator is infrared suppressed, whereas 
the ghost propagator is infrared enhanced. In the context of confinement, the results imply positivity violation 
(formalized in the Oehme-Zimmermann superconvergence relations [17]) and are in agreement with the Kugo-Ojima 
[18] and Gribov-Zwanziger [19-21] confinement scenarios. These results have been contested (see for example Ref. [22]) 
by appealing to the Slavnov-Taylor identity for the three-gluon vertex, applying the further assumption that certain 
dressing functions that compose the ghost-gluon vertex kernel remain infrared finite and with the conclusion that the 
ghost propagator is not infrared enhanced. The existence of such additional solutions has been verified in Ref. [23] (this 
type of solution was incidentally also observed in Ref. [24] within the canonical approach to Coulomb gauge Yang- 
Mills theory [25]): it was concluded however that the Slavnov-Taylor identities alone cannot discriminate between 
the two types of solution and that one needs explicit additional input (the authors arguing that the infrared enhanced 
ghost solution is preferred). 

From the above (admittedly historically incomplete discussion of) works, several lessons emerge. Firstly, the 
Slavnov-Taylor identities play a central role in nonperturbative Dyson-Schwinger studies since they provide informa- 
tion about the higher n-point functions that enter the Dyson-Schwinger equations. Second, in all the above studies, 
the Slavnov-Taylor identities relate the 4-vector contraction of a vertex to some combination of inverse propagators 
(proper two-point functions). This contraction means that only part of the vertex is constrained by the identity and 
it is not necessarily that part which enters the Dyson-Schwinger equations. One either assumes that the 'transverse' 
vertex (i.e., that part not constrained by the Slavnov-Taylor identity) can be neglected or considers more sophisti- 
cated input to complete the closure of the system such that the Dyson-Schwinger equations can be solved (e.g., the 
Curtis-Pennington vertex uses multiplicative renormalizability constraints [4]). Third, even if the Slavnov-Taylor 
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identities arc applied to close the Dyson-Schwinger equations, the solution (and in particular the infrared behavior) 
may not be uniquely specified and requires additional input, as demonstrated in Ref. [23]. These latter two points 
are emphatically not intended as a criticism of the Dyson-Schwinger approach, merely as reminders of some of the 
problems encountered (confinement is after all, not entirely trivial) and in this paper, we will see what Coulomb gauge 
has to say about the subject. 

The importancx; of Coulomb gauge in studying nonpcrturbative QCD was recognized early on, as was the inherent 
difficulty in technical calculation for such noncovariant gauges [26] . The significance of Coulomb gauge is based in the 
observation that in this gauge, the system reduces naturally to physical degrees of freedom (explicitly demonstrated 
in Ref. [21]). Given that the Slavnov Taylor identities arc the expression of charge conservation as applied to Green's 
functions, there exists a clear motivation to derive them in Coulomb gauge. In addition, there has recently been 
much technical progress in the Coulomb gauge functional formalism that provides the background to the study of 
the Dyson-Schwinger equations: the derivation of the Dyson-Schwinger equations themselves [27, 28], their one- 
loop perturbative analysis [28-30], the resolution of more formal aspects of the (incomplete) gauge-fixing and the 
existence of a conserved, vanishing total charge (and the absence of the infamous Coulomb gauge energy divergences) 
[31, 32]. With the Slavnov-Taylor identities, most of the components required for at least an initial study of the 
Dyson-Schwinger equations are in place. 

In this paper, we thus derive the Slavnov-Taylor identities of Coulomb gauge Yang-Mils theory in the second order 
functional formalism and discuss some of their immediate consequences. In Section 2, the basic formulation of the 
functional formalism and the Gauss-BRST invariance is presented. The Slavnov-Taylor identities for the two-point 
functions are derived in detail in Section 3. In Section 4, the identities for the vertex functions are derived and it 
is shown how these identities form closed sets such that the temporal Green's functions may be obtained as their 
solution. Section 5 is concerned with the possibility of extracting unambiguous information about Green's functions 
from the identities in the infrared. The paper concludes with a discussion of various aspects of the identities. For the 
convenience of the reader, selected relevant results from earlier works are listed in Appendix A. Lengthy configuration 
space expressions, important to the derivation of the Slavnov-Taylor identities but not to the narrative of the paper 
are relegated to Appendix B. 



2. FUNCTIONAL FORMALISM AND GAUSS-BRST INVARIANCE 



Let us begin by considering Yang Mills theory in the functional formalism. We will use the notation and conventions 
established in [27, 28] . We work in Minkowski space with, metric g^i/ = diag(l, —1). Roman subscripts {i,j, . . .) denote 
spatial indices and superscripts (o, b,...) denote color indices. We will often write configuration space coordinates 
{x,y, . . .) as subscripts where no confusion arises. 

The Yang-Mills action is defined as 



Sym = I d'^x 



4 /^^ 



(2.1) 



where the (antisymmetric) field strength tensor F is given in terms of the gauge field ^4^: 

F;, = d^Al - d.A- + gr'^AlAl. (2.2) 

In the above, the Z"*"^ are the structure constants of the SU{Nc) group whose generators obey [T°,T''] = if^^^T^. 
The Yang-Mills action is invariant under a local SU{Nc) gauge transform characterized by the parameter 6^: 

=exp{-^0»T'^} (2.3) 

such that for infinitesimal 0^, the gauge field transforms as 

A« ^ A'" =Al- -bfe^ (2.4) 

with the covariant derivative in the adjoint color representation given by 

D;^ = S''^d^+gr'''Al. (2.5) 

In terms of the temporal and spatial components, the above transform reads (we rewrite the temporal component 

as cr) 

9 

A«- A'" = > -I- -V^i" - P^^^e". (2.6) 
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Consider the functional integral 



J X>$ exp {iSy m 



} 



(2.7) 



where 2?$ denotes the functional integration measure for the collection of all fields. Since the action is invariant under 
gauge transformations, Z is divergent by virtue of the integration over the gauge group. To overcome this problem we 
use the Faddeev-Popov technique and introduce a gauge-fixing term along with an associated ghost term [33] . Using 
a Lagrange multiplier field, A", to implement the gauge- fixing, in Coulomb gauge (V-A" = 0) we can then write 



j V^expliSyM + iSfp}, Sfp = j 



(2.8) 



where and are the Grassmann- valued ghost fields. The new term in the action, Sfp, is invariant under the Gauss- 
Becchi-Rouet-Stora-Tyutin [Gauss-BRST] transform [21] whereby the infinitesimal, spacetime dependent gauge pa- 
rameter 6°^ is factorizod into two Grassmann- valued components: 0" = c%5\t, where 5\t is the time- dependent 
infinitesimal variation (not to be confused with the colored Lagrange multiplier field A"*). The Gauss-BRST transform 
is peculiar to Coulomb gauge — the time-dependent variation is allowed simply because the gauge-fixing does not 
involve any explicit time-derivatives. The variations of the new fields read: 



9 ^ 



5X1 = 0. 



By including a source term, the generating functional is given by 



Z[J] 



where 



/ 



j T)^ exp {iSym + iSfp + iSg} 



(2.9) 



(2.10) 



(2.11) 



Regarding the Gauss-BRST transform as a change of integration variables under which the generating functional is 
invariant, noting that the associated Jacobian factor is trivial [27] and only the source term varies, we deduce that 







-I 



5 [i6Xt 



■ exp {iSym + tSpp + ^Ss + iSSg} 



5\t=0 



j exp {iSym + tSpp + iSs} J d^xS{t — x 



- xa) X 



'K^xPxJ^x "T / Px^x^x ^"^ix^ix^x / "^ix-^ix^-x ^xHx q/ Hx^x^x 



(2.12) 



This equation is the starting point for deriving the Slavnov-Taylor identities. Notice the 5{t — xq) constraint, which 

arises because of the timc-dcpcndcnt variation SXt and is characteristic to the Gauss-BRST transform. It leads 
eventually to a nontrivial energy injection into the Slavnov-Taylor identities which is not present in the covariant 
gauge case. 

It is pertinent at this stage to discuss some nontrivial points associated with the above. Coulomb gauge is in fact 
not a complete gauge. Even after adding the gauge-fixing terms, the functional integral, Eq. (2.8), still contains 
zero-modes generated by purely temporal gauge transforms (and for that matter, global transforms too). Explicitly 
separating these temporal zero-modes within the Faddeev-Popov procedure, one can formally show (within the first 
order formalism, but since this is connected to the second order formalism used here via identities, the same conclusions 
apply) that the total color charge of the system is constrained to be conserved and vanishing [32]. The above 
expression, Eq. (2.12), is the dynamical statement of color charge conservation it shows how the external sources 
must be arranged in accordance with the gauge symmetry of the underlying theory. In the two differing contexts, 
the incompleteness of the gauge plays a central role: the formal isolation of the zero-modes gives the total charge 
constraint and the time-dependent variation (SXt) will give (later) the extra temporal scale. We will discuss this 
connection at the end. 

Note also that the presence of the temporal zero-modes precludes deriving a similar expression to Eq. (2.12) above 
by considering a purely time dependent (spatially independent) full gauge transform, i.e., 9{t,x) — > 0{t), since the 
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functional integration is ill-dcfincd. In such a case, the identity gives a resuhant Green's function (see later 
for details of the definition) with a dressing function that vanishes at zero momentum and this would hold even at 
tree-level which clearly contradicts the perturbative behavior for which the dressing function is unity (see Ref. [28] or 
Appendix A). The reason for this is that when considering the gauge- variant Green's functions, the integration over 
the gauge group associated with the zero-modes averages such quantities to zero. For the Gauss-BRST transform 
considered here, the spatial dependence of the ghost-field ensures that this problem is not encountered (since the 
gauge-fixing part of the action can be rewritten in terms of spatial derivatives of the ghost-field, they can always be 
implicitly defined as having no spatially constant component). 

So far, the generating functional, Z[J], generates all Green's functions, connected and disconnected. The generating 
functional of connected Green's functions is W[J] where 



We define the classical fields to be 



1 
z 



1 5Z 

Z SlJr, 



(2.13) 



(2.14) 



(we use the same notation for both the classical and quantum fields since no confusion will arise). The generating 
functional of proper Green's functions is the effective action, F, which is a functional of the classical fields and is 
defined through a Legendre transform of W: 



r[$] = w[j] - 



(2.15) 



In the above, we use a compact notation for the sources and fields: a generic field is denoted with source Ja 
and with the index a standing for all attributes of the field in question (including its type); further, summation over 
all discrete indices and integration over all continuous arguments is implicitly understood. We introduce a bracket 
notation for derivatives of W with respect to sources and of F with respect to classical fields (no confusion arises since 
the two sets of derivatives are never mixed): 



<iJa>- 



5W 

StJa ' 



6T 



(2.16) 



We can now rewrite Eq. (2.12) as 

= 



r 



Knowing the functional form of the ghost equation of motion [28] , 

<zc>= -Vlcl+gr'"^^i. [<^Jlvin> +AIc\ 

allows us to write 

V. 



(2.17) 



(2.18) 



j d'^x5{t - xq)^ <iA1^> Vixcl = j d'^x6{t - xq)^ 
= J d^x 6{t — Xo 



(-V2) 



<iAt> 



.(-V2) 

Inserting the above into Eq. (2.17) then gives 



(-V2) 



I}- 



(2.19) 



= J d^xS{t-Xo)i^-^ {d° <za» ct - r'"' <iu:> a\cl - i 



(-v^) 



1 



1 



1 



(2.20) 
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where tij(x) = Sij + ^ix^jx/i^^x) ^^'^ transverse projector in configuration space. 

Before considering the fuU content of the Slavnov- Taylor identities arising from Eq. (2.20), let us show that func- 
tional derivatives with respect to the Lagrange multiplier field, A", play no further role. The field equation of motion 
for A°, valid in the presence of sources reads [28] 

<zA>= -VixAl (2.21) 

such that the only non-zero functional derivative of <«A"> is 

<i4j^iA>= iS'-'^Vja^Siy - x). (2.22) 

All other functional derivatives of <iX^> vanish, even in the presence of sources. Taking the functional derivative of 
Eq. (2.20) with respect to lAf gives 



= J d'^x6{t-xo)^^S{z-x) <i-4> ~ 



(-V2) 



<iXiiA°'^> 



<ic't> 



r"'^ <iXliA^,> U^{S)A%cl - /""^ <iXliA1^> Uj{x) <tJli%> 



-f 



abc 



<ta"> 



^ <iplvrfx> -r'"' <%A1^> Uj{x)^ <iJlvrfx> +^/"'^ <^c> ^ <^r?^^>| .(2.23) 



(5iAf '^^ ' ''' 'StXf_ -"J-^-'^' '2-' SiX^ 

Now, taking partial functional derivatives of the Legendre transform, Eq. (2.15), one can show that 



<lJalJf3> <lJalJK><i^KiX^l^j><tJ.ylJf3> 

= 



(2.24) 



since all three-point proper functions involving functional derivatives with respect to zA^ vanish. Further using the 
result Eq. (2.22) we arrive at 



= J (fx5{t-XQ) 1^ 



6{z — x) 



^ ix^ iz 

Fvf) 



■5{x — z) 



<i4> -if""" [VizSix - z)] tij{x) [A%cl+ <iJ%vfg>\ I . 



(2.25) 

Knowing that Viz5{x — z) = —Vix5{x — z) and using integration by parts on the last term, the above expression is a 
trivial identity. Thus, even in the presence of sources, functional derivatives of Eq. (2.20) with respect to iXj, give rise 
to an identity from which no further information can be obtained and we can set the classical field A° = (except for 
within partial derivatives used in conjunction with the Legendre transform) and ignore the Lagrange multiplier field 
from now on. 

Equation (2.20) is Grassmann valued and to proceed, we must first take the functional derivative with respect to 
icj,. Since we will be taking further derivatives, fields/sources must be retained (with the exception of A as discussed 
previously). The subsequent Slavnov-Taylor identities are thus functional derivatives of 



= j d-^x S{t - a^o) <zc^«CT>] i4 - ^ [d° <ia^>] 6{ 



z — x) 



<i^x> +- 
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<i^icf> 



-r"'^ <iciiai> [<ipivni> -loiici] - r"^ <ia:> 



<ipWx> -i(tIS'^''S{z - x) 



<iJlir]l> -lA^J'^'^Siz-x) 



+ Micl> [<vtirfx> -ic^icl] - Ir'^ <zc> 



610^ 



<ttxtrj> -2S^''S{z - x)i 



(2.26) 



SLAVNOV-TAYLOR IDENTITIES: TWO-POINT FUNCTIONS 



The first of the Slavnov-Taylor identities is generated by taking the functional derivative of Eq. (2.26) with respect 
to la^ and setting sources to zero. The resulting equation reads 



= J d^x6{t-xo)l^-^ [d° <ialiai>\ 5{z + ^ 



<zc"zcf> 
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-r"" <ial,ial> ^ <iplvft> <^<^A«,> Ujix) <tJj',zC> j • (3-1) 

"^'-z ,7=0 "''^z ,7=0 J 

This equation is best expressed in momentum space. We define the momentum space two-point Green's functions via 
their respective Fourier transforms 

<iaZial>= J dke-^''-^^--^Tt{ko,k), (3.2) 

(similarly for Taa etc. and also for the propagators Waa etc.) where the arguments of the momentum space functions 
reflect their noncovariant nature, we impose translational invariance in the usual way and dk = d*fc/(27r)*. We also 
define 

I dke-'^-^^-^^tl%,{ka,k), (3.3) 
j dke-^''<^-^^tX-,,^{ko,k). (3.4) 



j=o 



j=o 



The latter two definitions and notations may at first sight appear somewhat artificial. In fact, as will be shortly 
justified, their exact form in terms of standard Green's functions is irrelevant since neither term will contribute. 
However, in the remainder of the paper, similar (but increasingly complicated) objects will occur so some explanation 
is in order. The Dyson-Schwinger equations arc derived as functional derivatives of the field equations of motion and 
this gives rise to the familiar loop expressions. For example (taken from Ref. [28]), 

j d^yd^zT^^}'^\y,z,x)<vrf,ir^l>^ j dujT^'^\uj-k,-uj,k)W£j{uj)V'^:^{.{uj,k-uj,-k)W^^^^ (3.5) 



gives the ghost loop term of the Dyson-Schwinger equation for F^^ (VF-c is the ghost propagator, V-^a is the spatial 

ghost-gluon vertex and Fi°^ is the tree-level counterpart). Importantly, all occurrences of objects such as <irf^ir]y> in 
the field equations of motion have an associated tree-level vertex - these terms originate directly from the interaction 
terms of the original action such that the functional derivatives have a clear meaning as loops. However, in the 
case of the (nonabelian) Slavnov-Taylor identities one has a different structure which arises from the Gauss-BRST 
transform (or indeed generally from the BRS transform) and is exemplified in Eq. (2.26). In this case, one has objects 
such as <ip\.ifi%> which must be functionally differentiated, but which are separated from any interaction type of 
factor. However, even without the interaction term, these functional derivatives do still have a partial meaning as loop 
expressions; it is simply that the tree-level vertex is missing (although the color factor and momentum conservation 
are present). We will denote such irregular expressions with a tilde (as in E above or F later) and since these terms are 
important to the arguments presented in this study, we will write out their explicit forms where necessary. As alluded 
to above, this type of pseudo-loop expression is common to the Slavnov Taylor identities of nonabelian theories and in 
linear covariant gauges, there is a familiar example, referred to as the ghost-gluon scattering-like kernel which appears 
in the identity for the three-gluon vertex [34]. 
Returning to Eq. (3.1), we can now write 

= / go ^fce-«'<'(*--'')-'*'-(— [-Ki{ko, k) - -^Tl\,(ko,k)T^i{qa + fco, fc) 

19 9 k^ 

-F- (fco, k)tl%{qo + fco, k) - K\,{ko, k)Uj{k)tf^^^{ko + qo, k)} . (3.6) 

In the above, notice how the 6{t — xq) factor characterizing the Gauss-BRST transform leads to the energy insertion 

go which breaks the energy flow through the various Green's functions and in particular how this insertion affects only 
the ghost functions. Since t, w and z are arbitrary, we can write down the momentum space Slavnov-Taylor identity: 



^Kiiko, k) = '-^KAiiko, kn^{qo + ko, fc) + F- (fco, fc)S«f,,(go + fco, k) +Tl%{ko, k)Uj{k)t'X%J 



Now consider the definition of I^Aj;ccj given by Eq. (3.4): under a parity transform, the vector source Jj*^ changes 
sign [27] and we see that 

^'X%ej{ko,-k) = -^l%ej{ko,k) (3.8) 
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from which we can infer that 

tij {k)t'X%^j{kQ + 90, k) ~ Uj {k)kj = (3.9) 

such that the last term of Eq. (3.7) vanishes. Because the energy scale, qo, is arbitrary, we can make a further 
translation qo —qo — 2fco and Eq. (3.7) becomes 

^Kiiko, k) = -^K\{ko, m^{-qo - ko, k) + r- (fco, k)tl%,{-qo - fco, k)- (3-10) 
9 9 k-' 

Knowing that the cr-field and p-source change sign under time-reversal, whereas the ghost- antighost pair is invariant 
[27], we have that 

K%c{-'lo-ko,k) = -E< (9o + fco,fc), 
ri^{-qo - ko,k) = T^J{qo + ko,k), (3.11) 

such that 

^Tliiko, k) = -^K%{ko, k)T%i{qo + ko, k) - TZiko, fc)%(9o + fco, k). (3.12) 
9 9 k^ 

Comparing this with the original Eq. (3.7) above, we find that 

t%^iqo + ko,k)=0. (3.13) 

The Slavnov-Taylor identity now reads 

^riiiko, k) = -^riXiko, k)Tt{q, + ko, k). (3.14) 

9 9 k-^ 

However, since qo is arbitrary, we can further state that the two-point ghost proper function Fee is independent of 
energy, i.e., 

TfM^ + ko,k)^T%i{k). (3.15) 

This was known to all orders in perturbation theory [27, 28] , but it is reassuring that the result is confirmed nonper- 
turbatively. We now have the final form of the first of the Slavnov-Taylor identities: 

^r^i(fco,fc) = -^VZM{ko,k)Vt{k). (3.16) 
9 9 k-' 

Repeating the above analysis but starting by functionally differentiating Eq. (2.26) with respect to lA^^ and setting 
sources to zero leads to the similar Slavnov-Taylor identity 

^Tt,iko,k) = -'^r^^Aki{ko,k)r^J{k). (3.17) 
9 9 k'' 

Now, since the energy, fco, is a scalar quantity we see immediately from Eqs. (3.16) and (3.17) that the two-point 
proper Green's functions involving functional derivatives with respect to the temporal gluon field, a, (these will be 
referred to as the temporal Green's functions) can be unambiguously expressed in terms of the ghost and (longitudinal) 
spatial gluon two-point proper Green's functions. (We will discuss the identities, Eqs. (3.16) and (3.17), and their 
general kinematical decompositions further in Section 5.) In a sense, up to the level of the two-point functions the 
field (T has been eliminated from the system or equivalently has been integrated out of the functional form of the 
action. The two identities, Eqs. (3.16) and (3.17), have been derived previously from the standard BRS invariance, 
using perturbative arguments to eliminate the S terms and were verified to one-loop order in perturbation theory 
[28]. The above derivation is however slightly superior since it makes reference only to symmetry arguments and 
does not rely on the implicit assumption that the all-orders resummation of perturbation theory is equivalent to 
the nonperturbative theory. Together, the two identities are the Coulomb gauge equivalent of the standard covariant 
gauge result that the longitudinal part of the gluon polarization is not dressed [1, 2]. This also highlights an important 
difference between Landau gauge and Coulomb gauge: in the former, the gluon polarization is transverse to the four- 
momentum; in the latter, the gluon polarization {Taa in our notation) is explicitly not transverse instead, the 
spatially longitudinal, ghost and temporal two-point proper functions are all related. This connection is reminiscent of 
the quartet mechanism in the Kugo-Ojima confinement scenario [18]. This is hardly surprising since the Kugo-Ojima 
confinement scenario is based on the identification of a well-defined conserved color charge and the Slavnov-Taylor 
identities here are the dynamical expression of charge conservation. 
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4. SLAVNOV-TAYLOR IDENTITIES: VERTEX FUNCTIONS 



In this section, we study further functional derivatives of Eq. (2.26) and show that the resultant Slavnov-Taylor 
identities eventually form a closed set from which the temporal Green's functions (i.e., those with at least one external 
cr-leg) can be unambiguously derived. This closed set turns out to be somewhat extended and so, for clarity of 
presentation, we separate the various classes of functional derivatives according to how many pairs of ghost /antighost 
functional derivatives are present. By explicitly separating the Grassmann- valued fields, we may then modify our 
notation for the generic field or source Jq. to be restricted to only the A oi a fields/sources and their associated 
attributes which allows us to compactify the formalism. 



A. No further ghost derivatives 

Let us begin by functionally differentiating Eq. (2.26) twice with respect to i^x(k)w ^^'^ *^t(0'« ^■'^^'^^ mentioned 
above, the indices A(fc) and t(Z) refer here to either the A-field with its associated spatial index (fc or I) or to the 
a-field with no associated spatial index. Setting sources to zero, the resulting equation reads: 







j <fx5{t- 



9" <iK(i)v''^mwi(^x> 



6{z -x) + 



1 



+- 



V,;, 



<»*A(fe)«,<> 



Kleins ta^> 



pabc 



t{1)v z 



abc 



J=0 



J=0 



r^'^S^xSiw - x)S{z - x) 



eabc 



J=0 



J=0 



r^%l6ArS{v - X)6iz - X) 



r%kSAxS{w - x)5{z - x) 



(4.1) 



In the above, we have made use of the results Eqs. (3.9) and (3.13) to eliminate such terms. Because the indices A and 
T refer to field types that commute, the above equation actually represents three separate equations. Let us define 



7=0 



be 



jb -ttc ^ 



(4.2) 



j=o 



(these expressions will be explained later). We further notice that by taking two functional derivatives of Eq. (2.18) 
and setting sources to zero then we have 



<ic:zcii<l>i^^^^>= -gr^'SArS.iVj.Siv - x)6{z - x) - V,.r^'^^^^.(,)(x, z, v). 



adf 



(4.3) 



This equation is of course the precursor to the Dyson-Schwinger equation for the temporal and spatial ghost-gluon 
vertices. Our identities can thus be written 



= J d'^x6{t-Xo)^-i 



^x <»^^(i).**A(fc)«,«^x> 



d{z -x) + 



ntcxm + gr^'S.kSAxSiw - x)5{z - X) 
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f :tcA(fc) z, w) + gr'''6„xS{w - x)6iz - x) 



(4.4) 



Let us now introduce our convention for the Fourier transform of the three-point functions (similarly for higher n-point 
functions): 



f{x, y,z) = J dki dk2 dkz (27r)''5(A:i + + kz)e 



ikx-x — ik2-y—tk^-z 



/(fci,fc2, Ais). 



(4.5) 



All momenta of the three-point functions are defined as incoming and for brevity we do not split the temporal and 
spatial sets of arguments. Returning to Eq. (4.4), after some manipulation we arrive at the three Slavnov-Taylor 
identities (fci + + fca =0): 



^3rflV)(fci,fc2,fc3) 



^ferf-wi(fci.fc2,fc3)r§,^(-fc3) 



■ \A{k)i 



(fc2) 



f;^L^(/)(^2 + go,fc3 



90 



'^A■-ccX^ik){kl + ao, fc3 - 50, ^2) + gP^'S^kSAy 



-niik) ik2) [f (fc2 + «0, fc3 - 90, fcl) + gr'^^. 



-r 



/a 



7(1) 



iki) 



cr;ccA(fc) 



(fci+go,fc3-go,fc2)+5r'''^aA 



(4.6) 



These identities can easily be verified at tree-level, using the Feynman rules presented in Rcf. [28] (and for completeness 
presented in Appendix A). They are the Coulomb gauge analogue of the familiar Slavnov-Taylor identity for the 
three-gluon vertex [34]. As in the previous section, the energy qq is injected into the ghost line of the various factors 
but one cannot cancel terms using symmetry properties anymore (as was the case for the S kernels) because of the 
presence of other energy arguments in the equations. Neglecting the F terms, we can see that each vertex involving a 
functional derivative with respect to the temporal a-field is fully determined (again because k^ is a scalar quantity) 
given the corresponding (spatially contracted) vertex involving the derivative with respect to the spatial A-field. 

Let us now discuss the form of the F factors. Prom the definition, Eq. (4.2), and taking functional derivatives in 
standard fashion, we have 



■pad/ 



.{x,z,v) = gf 



be 



6^ 



J=0 



(4.7) 



where we recall that having explicitly extracted the Grassmann- valued fields, the repeated indices refer here to only 
the a or A-fields. The coordinates m, s and r are implicitly integrated over. We further omit the possible spatial 
index (/) associated when r refers to the A-field for notational convenience. A similar expression exists for V a-;cct- In 
momentum space, we have (pi + P2 + P3 = 0) 



^adf 



.{puP2,P3) = gr'^J dkwll{k)Till{k,vz,-k-pz)w';:,{k+pz)w^^i{p^-k)^^^ 

-gr'" j dk Wl% {k)W^^ (pi - k)T'l^^^l {Pi-k,p2,k,p3) 

+gr'^ j dkW^Uk)Wi^{p, - k)r^^li{p, - k,p2 + k,p3)W^{-p2 - A:)rif^(-p2 - k,p2, k), 



(4.8) 

(we make the internal color indices explicit) again with similar expressions for the other F factors. Given that there is 
the implicit summation over e, jj and k (referring to the field types a and A), what this means is that the F considered 
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so far arc in general dependent not only on the non-gliost temporal three-point functions that arc explicit in Eq. (4.6) 
(thus showing that Eq. (4.6) actually forms a set of nonlinear integral equations) but on the further temporal Green's 
functions Tcca: ^ccaA and Tccaa (and of course all the other two-point functions, the purely spatial A and the ghost 
three- and four-point functions). 



B. One pair of further ghost derivatives 



Let us now consider the Slavnov-Taylor identities that arise if we functionally differentiate Eq. (2.26) with respect 
to at least icl and ic^. If we restrict ourselves to considering no further ghost /antighost functional derivatives, then 
we are able to set the corresponding fields/sources to zero whilst maintaining the rest without confusion. The resulting 
expression reads: 



= y d'^x6{t - a;o) 1^ [d° <icii4ial>] 6{w - ^) - ^ [d° <«4 



|0 ^,7;/,^e 



(-V2) 



.(-V2) 



<iAf^> 



<iK<n> +- 



(-V2) 



<iclic^iA^^> 



<t^ic'^> 



<^C^ ICyj IC^ ZC^ > 



5^ 



SiclSicljSicf 



5 



<ipli%> -i(jl5'^''5{z - x) 



<iJUr]l> -iA]j^'^5{w - x) 



+^r'" <4<> 



5^ 



^ <iJ%i%> -iA]j'^''5{z - x) 



5icf„5ici 



<irti%> -25'"'5^''5{z - x)5{w - x) 



r'"'<iAt,> Ujix) 



rf—ij—O 



(4.9) 



To derive the identity for the ghost-gluon vertex, we further set the remaining sources to zero (we will return to 
Eq. (4.9) later in this subsection to derive more identities) and again using the results given by Eqs. (3.9) and (3.13) 
we get 

= J d'^x6{t-Xo) [d° <icii4ial>\ 5{w - - ^ [9° <i4icltai>] 6{z - x) 



(-V2) 



<ic(,ic^iA^^> 



(-V2) 



<iclic^iA^^> 



+ ^ <ZC^C> fg.tee(^, Z, - r'^' <icilC> 5{z - x)5{w - x)^ 



(4.10) 



In the above, we have defined 

^ade 



V^,,{x,z,w) = gf 



abc 



Sicf,.Sici 



J=0 



= gf"'"' <iftir]y><iff^i-n^> [2 <zJ«zJe><zc^zcf2$«><«c^«c^z$e> - <«c,,i<icf?c^>] (4.11) 

(internal indices n and e referring only to A or cr-fields as in the previous subsection). Considering now the Fourier 
transform and using the same conventions as before, we arrive at the Slavnov-Taylor identity for the temporal ghost- 
gluon vertex 







/C3°r4t(fci,fc2,fc3) ^ ^!^r4''«i(fci,fc2,fc3)r--(-fc3) 



-{kl + 9o)r|!'(A;i, fcs - 90, fc2 + 3o) + ^r|!l(A:i, fcs - 90, ^2 + 5o)r^,^(-fc2) 



k'2 



+^r4"(fci)f^;tcc(Ai + 90, fc2, ^-3 - go) - gr'-ri^Ak) 



(4.12) 
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with 

^^{Pi,P2,Ps) = gr'^I dkw^jip, - k)w-:^{k) {2r4tbi - k,P2,ps + fc)r§::(fc,P3, -PS - k)w!:U-Ps - k) 

-^^:e'iPi-k,Ps,P2,k)}. (4.13) 

Recall that we are free to choose the energy injection, qo, at will. If we set + 90 = in Eq. (4.12), then (given 
that for such general spacelike momentum configurations, the Green's functions can have no singularities) the Tcca 
term of the middle line drops out and we have an unambiguous nonlinear integral equation for Tcca- Crucially, as 
far as the three-point functions involving the temporal cr-field are concerned, the equation forms a closed expression 
without approximation. One can immediately verify this identity at tree-level using the Fcynman rules of Ref. [28] 
(and summarized in Appendix A). A similar Slavnov-Taylor identity for the ghost-gluon vertex in Landau gauge does 
in fact exist [6, 12]. 

Returning to Eq. (4.9), we take the further functional derivative with respect to ^'^x{k)u ('^S'^i^i the field type A 
refers only to either the spatial A-field with index k or to the temporal a-field) and set sources to zero. Making 
use of the results Eqs. (3.9) and (3.13) to eliminate such terms and if we further use the appropriate functional 
derivatives of Eq. (2.18) as before, then we have Eq. (B.l) (such lengthy configuration space expressions are relegated 
to Appendix B). Equation (B.l) involves three new kernels: r<7;cccc and FqcccA are given by Eqs. (B.2) and (B.3), 
respectively; the third kernel, Ta.-ccc, has a similar expression to Eq. (B.2). There is only one possible new Green's 
function involving derivatives with respect to the o--field and this is Tzccca- Taking the Fourier transform of Eq. (B.l) 
leads to the following Slavnov-Taylor identity in momentum space (^ ka = 0): 







klTf^Zik){kiMMM) - =;r^iTAik)iiki,k2,ks,k4)m-k4) 

ih 



tk^i 



^fdga 



{kl + q,)Ti2' (fci, /C4 - 50, ^3, k2 + go) + -^^:iMk)iikiM - Qo, fca, k^ + Qonchh) 



^2 



+^-ccaikl,k2,ks + k4) 

+HcMikuk2,k3 + k4) 



r:r^Aw(9o-fc3 



A[ccXi(k) 



^fea 



^icAiiki,k4 - qo,k2 + ks + qo) 



(fci, ki - qo, k2 + k3 + qo) 



4r: 



'A:(fe)(fc3)f :%,(/c3 + go, k2, fc4 - go, fci) - rl\^k)iik^)^ZL-ciiks + 90, k2, k4 - qo, fci 



k4,k4-qo,k3)-gr^''S^x 
ks - k4, k4 - qo, k^) - gf^'^'^hi^AX 

,^{~k2~k^,k2,k^)-grH^x 
A;cc\i(k)(-'^2 - ^3, k2,k3) - gr<"^6kidAx 

adef 



adg 



o-:.ccX(k) 
^adg 



fag 



2 ccA(fc)^ 



(fci , k2 + k4, ks) T^cMo -k2- k4, k2,k4- qo) - 2gf 



ade 



-If 



fa 



(^i)f-:;?/cA(.)(fci 



qo,k2,k4 -qo,k3). 



(4.14) 



with the kernels (because of the proliferation of color indices, we resort to Greek superscripts) 

Cgc(pi'f'2,P3,P4) = gr'^J skw^:{k)wf](pi -k)x 

{ri':;(P4: -P4 - k, k)W^^{p4 + k) [ntiPi + k,P2, -k - P2 - P4)W:^{k +P2+ P4)q::(pi - k,P3, k+P2+ Pi) 



-r|S(P4 + k,p3,Pl +P2- k)W^;:{k-p,-p2)T^^:{pi - k,p2, k-p, -P2) 



+4^:ip4, k + P2,Pl +P3- k)Wl"{-p2 - k)T^i^{-p2 - k,P2, k)WI^{k 
^4^:iP4, k + P3,p, +P2- k)W-^;{-p3 - k)TB--{-p, - k,P3, k)W^^{k 
+r^3(P4,P3,-fc -P3 -P4,fc)W^r(fc + P3 +P4)r2^t(pi - k,p2,k+p3 
P2-p4,k)W^:{k+p2+P4)r2',^{pi - k,P3,k 



P2 



+ P2+ P4)r2'"(Pi - k,P3, k+P2 

+ P3+ P4)r2^t'(Pl " k,P2, k + P3 
+ P4) 
+ Pi) 



Pi) 
■Pi) 



cccc 
cccc 
ccv 



+ri7c (P4, k + P3,P1 - k,P2)Wl"{~p3 - k)T^:::i-p3 - k,P3, k) 
ip4,k+p2,Pl - k,P3)W^,"i-p2 - k)T^^:i-p2 - k,p2,k) 



+rir (P4, -P4 - k, k)w-^:{p4 + k)T;z^{p4 



k,P3,Pl - k,P2) - ^il^c^ij>4,P3,Pi - k,P2,k)'^ , 



(4.15) 
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r-:;?/ca(pi'f2,P3,P4) = gr'^'l dkwi^{k)w-:^{p^ - k)x 

{2r^^^j'(fc, -p, - k,p,)wl'{p4 + fc)rg2b4 + fc,P3,Pi +P2 - k)w::{k^Pi-P2)r2j^iPi - k,p2,k-p,-p2) 

-2^^^{k, -P4 - k,P4)wi%p, + k)nt{p4 + k,P2,pi +P3- k)w::{k - Pi - P3)r;::{pi - k,P3, k-p,- p^) 

-2r^"'(fc.P3, -P3-P4 - fc,P4)r2^tbl - k,P2, k+p3+P4)W^^{k+P3+p4) 
+ 2r^ct'(^'P2, -P2-P4 - k,pA)Tl2{pi - k,P3,k+p2+P4)W^^ik+p2+P4) 

+2r^-(A;,p3, -P3 - k)W^li{p3 + fc)r;i^«(p3 + k,p4,pi +P2 - k)WZ{k-pi -P2)Tlt:{pi - k,p2,k-pi-p2) 

-2T';^l{k, -P4 - k,p4)Wi^{p4 + fc)r^:J/(p4 + fc,P3,Pl - fc,P2) + r^:Jf/(fc,P3,Pl - k,P2,PA)} ■ (4.16) 

(A similar expression exists for the kernel Taicccc-) A-S before, since we may choose go at will and the energy (^4) is 
scalar, we now have two more expressions that are able to give T-caA and T-caa in terms of all the other previously 
considered Green's functions involving the cr- field except one: T-c-ca- Also as before, one can immediately verify this 
identity at tree- level using the Feynman rules of Ref. [28] (and repeated in Appendix A). 



C. Two pairs of further ghost derivatives 

In order to close the system of Slavnov-Taylor identities, we must find an equation for T-c-ca- This equation arises 
by functionally differentiating Eq. (2.26) with respect to ic^, icl, ic^ and ic^ and then setting sources to zero. Using 
the results Eqs. (3.9) and (3.13) and the appropriate functional derivatives of Eq. (2.18) we obtain Eq. (B.4) (the 
new kernel will be discussed below). This expression is cyclic symmetric in the three ghost derivatives id^, ic^ and 
ic^ and is antisymmetric in the two derivatives icl and ic^. In momentum space, the identity reads: 



= klT'^g^:{kuk2, ks, fc4, fc5) - -^r^jLt{k^,k2, ks, fc4, fc5)r--(-fc5) 

kc 



+{kl + 9o)4'Jca (fci, A:4, fca, ^5 - go, + qo) - ^r^clcAi(^i' ^=4, fe, fcs - «o, + go)r^f (-fe) 

+{kl + qo)T^J^fikuk, - qo, ka, k2, k^ + qo) - ^r^://^«.(fci, A^s - go, ^3, k2, k^ + qo)TI^{-k) 

+r|cCT (^3, ki, -fcs - A:4)f^^f '|,-(go + fcs + k^, k^ - qo, k2, ki) 

+r|cCT (^3, ^5 - go, go - ^3 - ^5)f ^fcttc(^3 + fcs, k2,k4, fcl) 

+r|f;(A:3, fe, -fc2 - ks)tl%^Jqo + k2 + fc3, ki, k5 - go, fcl) 

+HcAiik3,k4, -ks - fc4)f X|ccci(90 + fcs + fc4, fcs ~ go, fe, fcl) 
+r|cli(^3, fcs - go, go - fc3 - ^5)? ij^ccs(^3 + fcs, fc2, fc4, fcl) 
+r|f^,(A:3, fc2, -fc2 - h)TA%cc-AQO + fc2 + fcs, fc4, fcs " go, fcl) 

-r|^^(fci, A:4, -fcl - fc4)f^^|j'^_(go + fcl + fc4, fcs - go, fc2, fcs) 
-ri^;(fci, fc5 - go, go - fcl - fc5)f :f5_(fci + fcs, fc2, fc4, fcg) 

-r^f^(fci, fc2, -fcl - fc2)f "^|£-(9o + fcl + fc2, fc4, fcs - go, fcs) 
-rt,^^,(fci, fc4, -fcl - fc4)f ^;iL,(go + fcl + fc4, fcs - go, fc2, fcs) 
-r^lKfci, fcs - go, go - fcl - fcs )f^^^^^-^ (fcl + fcs, fcs, fc4, fcs) 

~^^cAi(kl,k2, -fcl - k2)TA--cL-iilO + fcl + fc2, fc4, fcs " go, fcs) 

+ l^-cc-cciki,k4, fcs, fc2 + fcs) [f|l,(go - fc2 - fcs, fcs - go, fc2) - 2gr^3] 
+ 1^-Jciki,k5 - go, fcs, go + fc2 + fc4) [f5;|l(-fc2 - fc4, fc2, fc4) - 2gr'^] 

+ 1^-Jl"{kl,k2,k3,k4 + fcs) [f^;fec(90 " fc4 " fcs, fc4, fcs " go) " 2ff/'"^'=] 
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+ ^r4"(fc3)f|t?cc(90 + h, k4, h - qo, fc2, fci) - ^Tl^ihn^SciQo + kuk4, fcs - 9o, ^2, fcs)- (4.17) 

Aside from the Green's function that we wish to calculate {Tcccca) i there is only one further unknown kernel: T-^-cccccj 
given in Eq. (B.5). Importantly, this kernel introduces no new Green's functions involving functional derivatives with 
respect to the cr-field and in momentum space, it reads: 

f|t/i.(Pl,P2,P3,P4,P5) = gr'^'l Sfc - k) |-ir^f/|f (/C,P4,P5,P3,P1 - k,P2) 

+ 2T(jf{p5 ,k+p3+P4,Pl- k,p2)W-^^' {-k-p3- P4)r^f<f {k,p4,-k-p3-p4,P3) 
+ 2r4tf (P5,P2,Pl - k,p^, fc + p3)VK:|(fc +P3)r^r (fc,P3, -k - P3) 

+ 2r4';-'^^(p5,P2,Pl +P3 - fc,fc+p4)r^r(fc'P4,-fc -p4)W^„t(fc +P4)W-^t(fc -Pi - P^W-cfs (Pl " k , P3, k - pi - p^) 
+2r|c''/f (P5,-A; -P2 -P5,k+p2)r!^^^{k,P2,-k-p2)r^cs{Pi - k,P3,k-pi - P3)Tl^^'^{k + P2 + P5,P4,Pl +P3-k) 
>^W^I^{k+p2)W^^{k~p,-p3)Wl%k+P2+P5) 

+2rl^^'{p5,k+p2 +Pa,Pi +P3- k)^^"(k,p2,-k - p2)r^^l{pi - k,p3,k-pi - P3)Yl%{-k - P2 - P4,P4, +P2) 

y^w:i{k + P2)w-f,{k - Pi - P3)w-:^{-k -P2- Pi) 

-2r4'}f (P5, -k-p3 -P5, k + P3)r^'^^^{k+P3+P5,P2,Pl - fc, P4)r^c^ (fc, P3, -fc -Pa) 
>^W:^{k+p3)W-^,%k+p3+P5) 

+2r4l;c'^(P5, k+p3+p4,pi - k,p2)rl^^^{k,p3, -k - P3)Vtf^{-k - P3 -P4,p4,k+P3) 
xM^„t(fc + P3)W-^^{-k -P3- PA) 

- '^Hczf {P5,k+P3+P4,Pi ~ k,p2) r^f^ ( A: , p4 , - fc - P4 ) r|^|^ ( - fc - P3 - P4 , P3 , A; + P4 ) 

xW^l^{k+P4)Wl'{-k-p3-pA) + C.p. {C''{P2),C'{P3),CS{P4))} (4.18) 

where wc utilize the cyclic symmetry (cyclic symmetric terms denoted by c.p. and this includes the first factor F-c-,-,, 
which is the origin of the factor 1/3). Because of the cyclic symmetry and the fact that qo is arbitrary, Eq. (4.17) can 
be solved for Tcccca- and we have finally managed to close the set of Slavnov-Taylor identities. The identity, Eq. (4.17), 
trivially has no tree-level form. 

We have thus shown that just as for the temporal two-point functions, all temporal Green's functions considered 
so far can be deduced (at least in principle) from their non-temporal counterparts as solutions to the Slavnov-Taylor 
identities. The temporal a-field has been effectively eliminated from the system, or integrated out of functional form of 
the action. Again, as for the two-point functions we see that the Slavnov-Taylor identities relate temporal, (spatially) 
longitudinal and ghost Green's functions in a manner reminiscent of the Kugo-Ojima quartet mechanism [18]. 

D. Further Slavnov-Taylor identities 

It has been shown so far that there exists a closed set of Slavnov-Taylor identities that includes the vertex (three- 
point proper) Green's functions. However, one may also consider four-point fimctions (e.g., TAAAa)- Because of their 
extended nature, we do not attempt to derive these identities completely; rather, we shall merely sketch their form in 
order to justify that the equations close just as before. We begin as previously with Eq. (2.26) and again, the sets of 
derivatives are distinguished by how many pairs of ghost / antighost functional derivatives are taken. Wc will highlight 
only those temporal Green's functions or kernels that have not previously appeared and leave those quantities that 
have already been derived as implicit. The sequence is as follows: 

1. Taking functional derivatives of Eq. (2.26) with respect to i^r and (as before, the ghost derivatives will 

be made explicit and the field types here refer only to either the A or cr- fields), one clearly obtains an identity 
for r^rAcr in terms of T^txa and a new kernel T^.^cXt- This kernel, following (as a further functional derivative 
of) Eq. (4.7) introduces the new quantity: TccktX- Starting with Taaaa, we can then sequentially build up to 
^auaa, given all the TccktX- 

2. Next wc take one pair of ghost functional derivatives (i.e., functionally differentiate Eq. (2.26) with respect to 
ic and ic) and then derivatives with respect to i^r and This gives us equations for the TccrXa- in terms of 
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the TcctXA and wc have two new kernels: ra-;-ccc\ and f c;cccAt- These kernels follow from Eqs. (B.2) and (B.3), 

respectively, and both introduce the new function: TcccctX- 

3. With two pairs of ghost functional derivatives of Eq. (2.26), plus one further derivative with respect to t^x we 
generate equations for TccccXa in terms of TccccXA, again with a new kernel: Tc^cccccX- This kernel follows from 
Eq. (B.5) and introduces a final Green's function: TccccccX- 

4. We finally take three pairs of ghost functional derivatives to get an equation for T-c-c-ca hi terms of TccccccA- 
However, this equation does not involve any further kernels and the set of equations terminates. 

Now, all these equations have the same characteristics as in the previous subsections: the energy is a scalar quantity 
and there is the energy injection scale qo, such that one has a set of unambiguous (albeit nonlinear and extremely 
long) equations from which the temporal Green's functions may be derived, given the set of spatial and ghost Green's 
functions as external input. 

We conjecture that in principle, the closed sets of Slavnov-Taylor identities in Coulomb gauge may be extended to 

include all higher n-point Green's functions and that they may be solved to give exact expressions for all the temporal 
Green's functions. As may be appreciated though, how to prove such a general statement is not clear. One may 
formulate some ideas based on the following observations: Firstly, the scalar nature of the energy and the energy 
injection scale {qq) will always be present in the Slavnov Taylor identities of Coulomb gaugc^ courtesy of the particular 
Gauss-BRST invariance of the theory (this certainly improves on the situation in linear covariant gauges where one 
has contractions of tensors such that transverse parts cannot be directly deduced) . Second, because ghost derivatives 
must come in pairs, each further pair reduces the number of possible non-ghost internal lines within the kernels - the 
external ghost legs must be connected by an internal ghost propagator or a vertex with these two ghost legs. This 
restricts the number of possible internal temporal Green's functions such that eventually the set of identities closes. 
One can see the emergence of a characteristic pattern for the ghost functions necessary to form the kernels in the 
steps above: i.e., (1.) r^cKrA ^ (2.) T^cZcrX ^ (3.) T'-^cZcccX ^ (4.) 0. 



Given that the Slavnov-Taylor identities derived in previous sections relate the various types of Green's functions 
in an extended manner, it is pertinent to ask whether these relationships reduce in specific circumstances such that 
unambiguous statements can be made abcnit the b(^havior of individual Green's functions, in particular the two-point 
functions. The motivation is clear: such information (particiilarly in the infrared region) may be useful for determining 
how the confinement mechanism manifests itself (asymptotic freedom already being perturbatively established in the 
ultraviolet region of Coulomb gauge [28]) and also provides for input in other nonperturbative studies. We argue here 
that unfortunately such information cannot be unambiguously extracted. 

In the noncovariant Coulomb gauge, we must first sort out the temporal (energy) and spatial (momentum) scales. 
The most obvious simplification of the Slavnov-Taylor identities concerns purely spacelike momenta where we know 
that if the postulate of Euclidicity is to hold (i.e., that the Wick rotation is valid), the Green's functions presumably 
can have no singularities. One may in principle also consider timelike configurations, but this is certainly beyond 
the scope of the present analysis. Thus, we must first set all energy scales to zero. We define the infrared region as 
the limit as one or more of the momenta vanishes and we approach the light-cone, where singularities (or nontrivial 
zeroes) may appear. 

Let us begin by considering the Slavnov-Taylor identities for the two-point functions: Eqs. (3.16) and (3.17). Using 
the general decompositions of Ref. [28] (presented also in Appendix A) in terms of the (scalar) dressing functions, we 
get for general kinematical configurations: 



5. 



SIMPLIFICATIONS OF THE IDENTITIES 



-ikok'^r^Aklk^) - -ifcoPrA.(fc2,P)ree(P) 

ikQki~Tjif^[kQ,k^) = tkQkkTAA{kQ,k'^)Tcc{k^). 



(5.1) 



Since both equations have common kinematical prefactors, we can then write 



(5.2) 



This holds for any kinematical configuration and we conclude that the Slavnov Taylor identities for the two-point 
functions alone can give no information about the value of the Green's functions, only the relationship between them. 
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Slightly less trivial is the Slavnov-Taylor identity for the ghost-gluon vertex, Eq. (4.12). Concentrating on spacclike 
momenta as discussed above, we see that neither of the Fcca- vertices will contribute when the energy is set to zero. 
To study the equation, let us firstly decompose the ghost-gluon vertex into tree- level and dressed parts as follows: 

ncMiPl,P2,P3) = -igr'^'PU - ^g.r'"'Pl,rA-ccA,^{Pl,P2,P^). (5.3) 

Let us explain this decomposition. Since we are dealing with a three-point function, the color factor (Z"*"^) and 
coupling [g) will always be common and can be extracted. The appearance of V stems from the momentum space 
form of Eq. (4.3) and in fact, given the definition Eq. (4.8), the above equation, Eq. (5.3), is the Dyson-Schwinger 
equation for the spatial ghost-gluon vertex. As will be seen shortly, the contraction of pij with r^.^cAji is important. 
Lastly, for the F^iccc kernel, we can only extract the color and coupling factors and we write 

f|'.ee = 9r''r-,-,,c (5.4) 

At zero energy, Eq. (4.12) in terms of dressing functions is thus 

= -ki-k^Tccikl) - ki-k2Vzc{kl) - klVcdkl) - kijk3irA;ccAji{ki,k2, k3)rcc{kl) - kijk2irA;ccAji{ki, h, ^2)r5c(fc2) 

+ ^fc?rec(fc?)f e;ecc(fcl, fc2, h). (5.5) 

This equation is directly analogous to the Landau gauge Slavnov-Taylor identity for the ghost-gluon vertex which 
was studied under the truncation Pqccc ~ [-^^]* Quite generally, it states that the kernels (i^A;ccAji 

and Fcjccc) and 

the two-point ghost dressing function Tec are nontrivially related. 

Let us now discuss the above identity, Eq. (5.5), in the infrared to assess whether it simplifies further. Because of 
the symmetry, there are two limits of interest: fci ^ and k2 0. For the infrared limit ki 0, we see that the 
entire equation has an overall factor ~ |A:i| (in fact, one can show that there is the overall factor |A:ip). Therefore, 
in the limit, no single Green's function is isolated from which to determine a value for this kinematical configuration. 
Further, the Green's functions or their combination may even be singular. For the limit k2 — > 0, we recall that in 
Coulomb gauge, ghost vertex dressing functions vanish as the "in-ghost" momentum vanishes [27], i.e., 

^A;ccAji{kl,k2,k3) OC ° \k2\, 

Tc-ccciki,k2,k3) '"oT" \k2\- (5.6) 

This is, of course, exactly the same as in Landau gauge [2]. Using the momentum conservation, fca = ~ki — k2, and 
we see that Eq. (5.5) has the overall factor \k2\ as fc2 0. This leads to the conclusion that no information about an 
individual Green's function may be extracted in this limit either. We conclude that on its own, the Slavnov-Taylor 
identity for the ghost-gluon vertex, Eq. (4.12), does not yield unambiguous information about the Green's functions 
without further knowledge. 

For the Slavnov-Taylor identities given by Eq. (4.6), there is one special case where simplification occurs and this 
again results in a situation exactly analogous to Landau gauge. Considering the Faao- identity (i.e., setting the 
external indices A and r to be referring to A-fields), one has explicitly 

k°sr{:Lkiki,k2,ks) = i'^Ti^z,,{ki,k2,ksn^{-h) 



k. 



-^TaM [r''A%cAiiik2 + qo, ks - qo, fci) + gr'^^du 
-^AAuiki) [n%Aik{ki + qo, ks - qo, fcs) + gr'^'Sik 

~rZkik2)K%,Aiik2 + qo, ks - qo, k,) 

-'^ti{ki)K%Ak{ki + qo, ks - qo,k2). (5.7) 

Further setting all energy scales to zero and using the general decompositions {T'^a — —Wf°'^'^^3A, for the rest, see 
above or Appendix A) one has 

k3i^3Alki{ki,k2, k3)rcc{kl) = kltkl{k2)TAA{kl) + kjtki{k2)TAAik2)'rA-;ccAzl{k2,k3, ki) 

-kltikiki)TAA{kl) - kltii{ki)TAA{kl)tA;ccAik{ki, fcs, k2). (5.8) 
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Taking the contraction of the above equation with k2k, one has thus 

k2khi^3Aiki{ki,k2,k3) = -kl ^^^j-^^\ ij{ki) fc2j + fefef ^;ccAi;o(4, 4, ^2) • (5-9) 

Now, the overall Bose-symmetry of the three-gluon vertex means that after extracting the color factor (/"'"^), the 

dressing function is antisymmetric under exchange of any two legs and so, by interchanging ^ k^i one can 
eliminate the three-gluon vertex to obtain an expression involving only the two-point ghost dressing function (r^c) 
and the kernel TA^ccAjk which, after canceling out the overall factors reads: 







tljiki) |fc2jrcc(fc|) + ksjT-^cikl) + Tccikl)k2ki^A^ccAjk{kl,k3, fe) + Fcc{kl)k3ktAiccAjk{kl,k2, 4)| . (5.10) 



In Landau gauge, this identity is well-known [35, 36] and was considered in Rcf. [22] under reasonable assumptions to 
be indicative of an infrared finite ghost dressing function (F-c). In the infrared limit A;2 — > (equivalently for fca — *■ 0), 
we see that all terms have the prefactor \k2\ since k3jtij{ki) = —k2jtij{ki) and for the kernel we have the general 
infrared result Eq. (5.6). Without further assumption, the following (conservative) statement is true of Eq. (5.10): 
if the particular contraction above of the ghost-gluon vertex kernel {J^A-;ccAjk) is vanishing, then the ghost dressing 
function is constant. However, notice that the 'particular contraction' above is not the same as that appearing in 
either the Slavnov Taylor identity for the ghost-gluon vertex, Eq. (5.5) or the Dyson Schwingcr equations from which 
Tec can be obtained. Thus we conclude, as previously, that the above component, Eq. (5.10), of the Slavnov-Taylor 
identity, Eq. (4.6) does not yield unambiguous information about particular Green's functions without further input. 

The special symmetric contraction case (above) of the three-gluon vertex Slavnov-Taylor identity aside, in order 
to say anything about the two-point functions from the identities Eq. (4.6), one must know something about at least 
one of the vertices Vr\a or Tr\A- Such information is not available from general considerations and so, one cannot 
make any simple statement. Further, the rest of the Slavnov-Taylor identities, equations (4.14) and (4.17), clearly 
involve too many higher n-point functions to have any hope of simplification. 

From the above discussion, it seems clear that no information about the behavior of the two-point functions (or 
vertices) can be obtained from the Slavnov-Taylor identities alone. As is obvious from the simple Slavnov-Taylor 
identities for the two-point functions summarized by Eq. (5.2), one can only determine the relationship between 
the various Green's functions and this must be true for any kincmatical configuration. Our only unambiguous 
'outside' knowledge about the Green's functions was the peculiar infrared behavior of the ghost vertices, the absence 
of singularities for spacelike momenta and the Bose-symmetry of the three-gluon vertex (which is why the ghost-gluon 
vertex identity, Eq. (4.12), and the identity Eq. (5.10) were of particular interest), but even then, the dimensionality 
or complexity of the Slavnov-Taylor identities denied concrete conclusions. 



6. SUMMARY, DISCUSSION AND CONCLUSIONS 



In this paper, the Slavnov-Taylor identities for Coulomb gauge Yang-Mills theory within the second order formalism 
have been derived. The starting point was the Gauss- BRST invariance of the action [21], characterized by the 
time- dependent BRS variation and peculiar to Coulomb gauge. It was found that for the two-point and vertex 
Green's functions, the Slavnov-Taylor identities form closed sets from which the temporal Green's functions can be 
unambiguously derived given the relevant spatial gluon and ghost Green's functions as input. The extension of this 
to higher order Green's functions was also discussed. Special cases of the identities were studied and it was seen that 
there is no simplification such that information about a particular Green's function could be extracted. 

It is worthwhile at this stage to discuss the Coulomb gauge Slavnov-Taylor identities in their wider context since 
the connection between several different themes becomes apparent. The first of these themes centers around Gauss' 
law. In classical electrodynamics. Gauss' law allows one to determine the temporal component of the gauge field (i.e., 
cr) directly from the physical charge distribution. In the canonical (Hamiltonian-based) formulation of Yang Mills 
theory. Gauss' law is applied as an operator identity to define the physical state space and can be explicitly resolved 
in Coulomb gauge, yielding the so-called Coulomb term which comprises the confining properties of the theory [25]. 
In the functional formalism. Gauss' law appears after the elimination of the temporal (tr) field from the action, either 
directly as in the case of the first order formalism [21, 32] or as here as the elimination from the effective action via 
the Slavnov-Taylor identities. The resolution of Gauss' law and its connection to the physical state space is certainly 
less apparent in the functional formalism since Green's functions are not directly related to physical obscrvablcs but 
this is compensated for by the observation that the cancellation of the gauge-dependent degrees of freedom (i.e., 
temporal, longitudinal and ghost) is manifest — in the first order formalism one can reduce the functional integral 
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to transverse spatial ghion degrees of freedom and the Slavnov-Taylor identities here express this expHcitly for the 
Green's functions in local fashion. However, the connection between the physical state space and the Green's functions 
is understood conceptually within the framework of the Kugo-Ojima confinement scenario [18]: by postulating a well- 
defined BRS charge and physical state space, the cancellation of the gauge-dependent degrees of freedom followed 
and in Landau gauge led to the prediction that the ghost propagator is infrared enhanced. In the Coulomb gauge 
functional formalism, we know explicitly that there exists at least a total charge that is conserved and vanishing 
[32] which partly confirms the Kugo-Ojima postulate. This total charge arises from considering the temporal zero- 
modes inherent to Coulomb gauge and leads to the formal demonstration of the cancellation of the gauge dependent 
degrees of freedom. In this study, we have shown that the Slavnov-Taylor identities stemming from the temporally 
nontrivial Gauss-BRST transform supply this cancellation in local fashion. In this respect, the temporal features of 
the Coulomb gauge functional formalism can be seen to supply a link between the physical charge and states to the 
Green's functions of the underlying theory. 

The temporally nontrivial nature of Coulomb gauge is manifested in two ways. On the one hand, resolving the 
temporal zero-modes leads (as mentioned above) to the vanishing and conserved total charge [32]; on the other 
hand, the extra temporal degree of freedom in the Gauss-BRST transform leads to the closure of the Slavnov-Taylor 
identities. The temporal zero-modes of the Faddeev-Popov operator lead us to consider the connection to the Gribov- 
Zwanziger picture of confinement [19-21]. In this scenario, it is recognized that the zero-modes (which are induced 
by incomplete gauge-fixing) should be separated from the functional integration and via stochastic quantization and 
entropy arguments, the authors were able to show that the resulting spatial transverse gluon propagator would be 
suppressed in the infrared (and thus drops out from the physical spectrum) whereas the temporal propagator provided 
for a long-range confining force. Just as with the Kugo-Ojima scenario, the ghost propagator in Landau gauge would 
also be infrared enhanced. In the Coulomb gauge functional formalism insofar as the temporal zero-modes are 
concerned, one can see the parallels: their resolution leads to a total charge and the cancellation of gauge degrees 
of freedom, providing an explicit demonstration of selected features of both the Gribov-Zwanziger and Kugo-Ojima 
confinement scenarios. However, one part is evidently missing — a prediction for the infrared behavior of the ghost 
(or something similar). Resolving the temporal zero-modes in Coulomb gauge restricted the functional integral to 
field configurations such that the total charge is conserved and vanishing, whereas the same temporal feature gave rise 
to Slavnov-Taylor identities that form closed sets but explicitly no information about a particular Green's function. 
This underlies the quite general feature of functional techniques, namely, that in order to talk about Green's functions 
one must take functional derivatives and in this respect, the Slavnov-Taylor identities (and for that matter, the 
Dyson Schwinger equations too) represent functional differential equations. Their solution is known only up to some 
'constant' of integration (for explicit examples of this, see Refs. [23, 24, 37] and references therein). In general this is 
not obvious since one necessarily has an infinite tower of equations to consider, but in the case of the Coulomb gauge 
Slavnov-Taylor identities derived here, the closure allows us to see clearly that this is precisely the case — functional 
techniques lead to relaMons between Green's functions. 

Pragmatically, there are two natural avenues which to explore further. The first is to search for the 'missing' 
prediction about the value of a specific Green's function (most likely the infrared behavior of the ghost) in order to 
complete the connection between the temporal aspects of Coulomb gauge and the two confinement scenarios: Kugo- 
Ojima and Gribov-Zwanziger. The second direction to take is to use the Slavnov Taylor identities here in order to 
construct a charge-conserving truncation scheme from which to solve the Dyson-Schwinger equations of Coulomb 
gauge Yang-Mills theory, allowing for the eventuality that the solution may only be determined up to some external 
'boundary condition' in analogy to the solution of standard differential equations. Both these directions are being 
pursued. 

On a final note, one further consideration for the Slavnov-Taylor identities derived in this paper is their verification 

to one-loop order in perturbation theory (the tree- level forms arc trivial). This has in fact been done. However, as 
can be appreciated from the lengthy expressions, such a technical exercise is not siiitable for inclusion in the present 
paper. The verification of the one-loop identities involves three components: the one-loop expansion of the vertex 
Dyson-Schwinger equations (and the kernels presented in the text), the use of the inherent translational invariance 
of the loop integrals and identities for the color factors. Importantly, none of the loop integrals need be explicitly 
evaluated. 



Acknowledgments 

This work has been supported by the Deutsche Forschungsgemeinschaft (DFG) under contracts no. DFG-Re856/6- 
2,3. 



18 



W Aj a__ _r Aj a 

A Uj{k)jj^^^DAA A -l{k^-k^)tij{k)rAA-tk^^TAA tk°kiTAa 

a ^D„„ a ik°kjTAa -ik'^T^^ 

A |i ^D„, A -ki 



TABLE I: General form of propagators [left] and two-point proper functions [right] (without color factors) in momentum space. 
All dressing functions are functions of ko and Jp . 



APPENDIX A: FEYNMAN RULES AND DECOMPOSITIONS 



For completeness, we present here selected results from Ref. [28] concerning the tree-level forms and general decom- 
positions of various Green's functions. The general decompositions of the non-ghost propagators and proper two-point 
functions are given in Table I (without the common color factor 6'^^). For the ghosts, we have 

wi'ik) = -s^^'d-uH rf,{k) = rtPr,,(P). (a.i) 

At tree-level, all dressing functions are unity. The trcc-lcvcl vertices are given by (all momenta are defined as incoming) 

C;%(Pa,P6,Pc) = l9r''SiM-P% 
r'-°tf{Pa,Pb,Pc) = -l9r'%Pa-Pc)j, 

^'ilijkiPa,Pb,Pc) = -igf^" [5ij{Pa - Pb)k + 5jk{Pb - Pc)i + hi{Pc - Pa)j] , 

^fJltiP^'PcPA) = -igr'^Pci. (A.2) 



APPENDIX B: LIST OF LENGTHY FORMULA 



In this appendix, we list lengthy (configuration space) formula, necessary for the derivation of the Slavnov-Taylor 
identities, that occur in the text. 



S{w — x) — I 

V. 



6{z — x) 



<iciic'iial> [f;:';|^^(,) {x, w, u) - gr^'S„x5{u - x)S{w - x) 
f :t A(fe) ix,z,u)- gr^XxS{u -x)S{z-x) 



+ <icii4iA1^> 



<tc{,ic^iA1,j,> 



K%cXi(k)(^' - ar'SkiSAxSiu - x)6iw - x) 

'^A%x^ik)ix, z, u) - g.r%iSAxS{u - x)S{z - x) 



1 



> 



^^,{x, z, w) - 2gr'''^S{z - x)S{w - x) 



<^Kik)u<> ^"J^-l-ci^^w)- <^Kik)u^^t> r'X%L{x,z,w,v) + \ <^4^c> ^%l,^^^{x, z,w,u)^ 



(B.l) 



^3 

r'"^f^-(x,z,w,v) = gf^" — J -<ip\irt> 



7=0 
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{<lcllCfj_l^i,><tri^ir]s><lJalJK> [<lCelc'^l^a><lC-ylC^l^K> — <tCstC^t^ aXiC-ylcfl^ k>] 

+ <ic{,iCfj_i^s><tri^iria><iJetJK> [<^Ca^c^^$l.><^c^^f•5;,?.<f>K> - <^Ca^c^^^'^/><^c-y^cf^$K>] 



+ <icliCKiCytc'!,><ir]^ir]a><iCaic'^i^,y> - <icliCK,ic^ic'^><vr]^iria><iCaic'^i^„> 

+ <iclic^ii^,y><vf]i^ir]K><tCKic'^iCjic'^> — <ic{,icljic^ic'^i^„>} . (B-2) 



abc 



J=0 



= gr^'' <triir]^><irf^ir]^> x 

-2 <iJsi^^> [<^Cy^c^^$£^$^^><^c-y^cf^$„> - <zc^zcfz$£Z<i>^^><zCyZc^i$K>] 

(In the above expression, Eq. (B.3), we omit the possible index (A;) when the field type A refers to the A-field for 
notational clarity.) 



= j d^x 6{t - xq) I [id^ Kic^iclici 



+ <^c^'^c^^c{^c^^c^f >] 6{u - x) 
+ <^c^Oc{^c>f>] 5{z - x) 



ic'^ia^>] 5{w — X 

V. 



+ <^4^4^a> f»g^(a;, ^, r)- <i4iclia^> tlg-{x, u, z, v) 
+ <^cl^<^<> z,w,r)- <i4iciia^> fl%{-{x,z,w,v) 

+ <iciiciiA1,> f%i^^-^(a;,«;,u,r)- <^c^cf^A^,> t''^$U{x,w,u,v) 
+ <tcliclzAt,> r%^,,-Jx,u,z,r)- <tc^,icltAt> f%{-^{x,u,z,v) 
+ <iciicitA1,> f''^g,-,i{x,z,w,r)- <ic'^iciiAl> f'X%i-Jx,z,w,v) 



r|;lc(^' - '^gr'Siw - x)5{u - x) 

rrccc(^> ^) - 2fl/'^^'^<5(« - x)5{z - x) 



f^%,{x, z, w) - 2gr^-d{z - x)S{w - x) 



-- <c^c> n^::dx,z,w,u,r) - - <rc-zcl> T-™ (cr, ^) 



(B.4) 



> 



5/"^'= <vff^ir]ij,><iff^tr]^> {- <tCi^iclcltcl,ic^i4> 

+2<if]^tr]s> [<icliCi^iCuic'l><ic^icliCsicl,> +c.p. (cf,c^,4)] 

+2<iJaij0> [<mli4K^icli<^s><ic^,ic';^,i^,,> +c.p. (cf,c^,cS)] 

+2 <lJalJi3><lJ-ytJs> [<iclic'^^i^^i^fj><ic^ic^i^a><Ki,ic'^i^s> +c.p. (cf,c^,c£)] 

+2 <ir]^iris><iJaiJf!><iJ'yiJs><icliCKi^fj> [<iCfj^ic'!,i^a><tc,yic'^i^s><ic^ic^i^-y> +c.p. (cf,c^,c£)] 

+2 <vr]^irjs><iJaiJp><iJjiJs><iclic^i^^> [<iCf,tc'^^i^a><iCi,ic'^t^s><iCsicli^p> +c.p. (c^,c^,0] 

-2 <Vq^irie><lJalJ0><cilCf,l^l3> [<lCelcfic„lcl><lC^lcl^l<^a> +C.p. (cf,C^,C^^)] 
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+2 <irj^irii;><iJaiJfj> X 

[<ic.liCi^iCvic'l> (<^c^^c^^^>a><^C£^c^^^>/3> - <ic^icli^a><iCsicl,%^f}>) + c.p. (cf , c^, c^)] } . (B.5) 

Note that in this expression, terms cyclic symmetric in the three ghost derivatives icf, ic^ and ic^ are denoted 
c.p. (c_j,, c^, c^). 
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